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Phase screens above a telescope pupil represent the variation of the phase of the electromagnetic 
field induced by atmospheric turbulence. Instances drawn from such statistics are represented 
by a vector of random phase amplitudes which are coefficients of a linear superposition of two- 
dimensional basis functions across the pupil. This work shortly reviews Fried's analysis of this 
modal decomposition for the case of Kolmogorov statistics of the phase covariance as a function of 
separation in the pupil plane. 

We focus on the numerical synthesis of phase screens. The statistically independent modes are 
transformed into the eigen-modes of a gradient matrix as time-dependence is introduced such that 
on short time scales the instances of the phase screens are rigidly shifted into a direction imposed 
by some wind velocity — known as the Taylor frozen screen approximation. This simple technique 
factorizes spatial and temporal variables and aims at binding the time dependence of the phase 
screens to the few expansion coefficients of the basis functions that obey a stochastic time-dependent 
' differential equation. 
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dl' I. OVERVIEW 
O 

^-j ' A. Phase Screen Snapshots 

^ ' A phase screen summarizes the phase delays in sub-apertures of a telescope of some diameter D on lateral spa- 
tial scales that depend on the inhomogeneity of the refractive index along the path through the few kilometers of 

'nI" • atmosphere in ground-based observations |24| . 

^ , The intent to remove the time-dependent sparkling of star light in (almost all) applications of astronomy leads 

' to adaptive optics and is increasingly important as telescope diameters increase and/or anisoplanatism puts limits 

J — on coherency in wide-angle observations |l5l. |40|. [s^ . The equivalent variables in long baseline interferometry spawn 

_^ interest in the definition of outer scales. 



o 



The impact on imaging would be determined by the degrees of freedom and stiffness of response in any subsequent 
adaptive optics correction, and by additional weighting with spatial filter functions as found in some fiber conductive 
foUowup systems 26, 33, 45, 50]. 
. This work is a contribution to the numerical simulation of phases (^(r, t) in some pupil plane taking a two-dimensional 
L| ' spatial vector r and a time t as arguments. We follow Fried's analysis of the build-up of the phase at some snapshot 
. ^ [ in time if the covariance (or structure function) follows a Kolmogorov power-law as a function of the distance between 
' two points in |r| < D/2 [3, '20'|. Assuming that the number and size of the speckles is founded on Gaussian statistics 
'\ after a sufficiently long path through the atmosphere, the phase can be synthesized by linear superposition of two- 
. . . 1 dimensional basis functions multiplied by a vector of random variables with individual Gaussian statistics. The 
statistically independent basis functions are constructed as eigenfunctions of the generic Karhunen-Loeve (KL) integral 
equations which incorporate isotropy and strength of the covariance in the integral kernel and the sampling limits of 
the statistics (here: the circular pupil) in the region of integration. 

B. Taylor Movies 

With some assumption of ergodicity, independent snapshots of phase screens are generated by throwing dice for each 
new set of expansion coefficients with a random number generator. On short time scales, however, the speckle boiling 
time (coherence time) would be overtaken by the assumption that the phase screen would move — in the fixed pupil 
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coordinates — rigidly and predictably as drawn by some wind velocity v |38l . |43| | . The expansion coefficients follow a 
stochastic Ornstein-Uhlenbeck-type differential equation which ensures deterministic (first-order smooth) alignment 
of the pictures of the Taylor movie, but allows transformation of each expansion coefficient within the amplitudes set 
by the eigenvalue of its KL mode. 

The alternatives to this ansatz are the generation of static phase screens much larger than the pupil diameter, 
from which a moving disk is cut off to define the region of interest 0, [3, [2l|, [2^, [H, 0, 53, IZH . The principle of 
embedding the pupil into a larger virtual phase screen will also be the guideline to define on which tiniescale how 
much randomness will be mixed into the expansion coefficients. The main difference is the efficiency of maintaining 
either a small set of time-dependent expansion coefficients with the option to generate movies of arbitrary duration 
(here), or maintaining a spatial stripe of the sky in computer memory which has the width of the pupil but a much 
larger length determined by the product of movie duration and wind speed. 

A summary of the well-known algebra of KL eigenmodes is provided in section |TT1 The differential equation and 
dynamic matrix of the motion implied by the Taylor frozen screen assumption is written down and discussed in section 
IIIH and notes of expanding this into a stochastic equation follow in section [iVl Some of the mathematical aspects are 
outsourced to the appendix. 



II. KL EIGENFUNCTIONS OF KOLMOGOROV COVARIANCE 



A. Separation in Radial and Azimuthal Functions 

A summary of the established theory and mathematics of phase patterns [l^, H^, [5l| is given in this section, 
accompanied by errata in Appendix [D] for the key papers. The phase of the electromagnetic field in the pupil plane 
is decomposed into basis functions F{r) and ffuctuating expansion coefficients a 

V.(r)=^a,i^,(r), \r\ < D/2. (1) 

j 

If the a are supposed to build a vector of independent scalar variables which vary from one phase screen to another, 
the covariance statistics between two pickup positions of the phase translates into a requirement of the basis functions, 
which are forced to become KL eigenfunctions of a linear integral operator. If the two-dimensional region of sampling 
the phase statistics matches the symmetry in the covariance matrix, a separation ansatz within the Fi is fruitful. 
In particular, if the domain of the pupil is circular and the covariance isotropic (only depending on the separation 
|r — r'l), the eigenfunctions can be separated in azimuthal and radial terms, 

Fj{r) ^ K!^'^Hx)Mg{0), 0<x, O<0<2tt. (2) 

If the phase structure function V obeys a power law 

I? = 2c^(|r-r'|/ro)'+^ 7 = 2/3, (3) 

one can essentially reduce the KL integral equation to a universal form by decoupling the scaling parameter D/ro 
from the Fj and moving it into the coefficients aj . The scale factor is [l^, |40| 
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1 + 7 



(l+7)/2 

2[(24/5)r(6/5)]^/^ « 6.883877. (4) 



The Fried parameter tq is essentially proportional to the 6/5-th power of the observing wavelength [20j if the wave- 
length dependence of the structure constant of the refractive index remains negligible. I shall work with basis functions 
which are normalized over the scaled distance x = 21a;|r|/£) from the pupil center, 

rt^ 1-2-K 

F^{Y)(fr= / a:if^«) {x)dx = / Af2(6l)d6l = 1. (5) 

D Jo Jo 

The constant is set to 1/2 if radial distances are measured in units of the pupil diameter D, and set to 1 if radial 
distances are measured in units of the pupil radius D/2. The only purpose of this notation is to keep track of both 
choices that have been in use in the literature. The azimuthal basis functions 

M,(0) = ./I^x|'^°f „^J' (6) 
' y 271 1 sm{q0), ^ ^ 
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will be labeled with positive q for the cosine type and negative q for the sine type, akin to the nomenclature for the 
two symmetries of Zernike polynomials [13] ■ The Neumann factor e is defined as in the literature on Bessel Functions, 

, 2, \q\ > 1. 

The radial eigenmodes for each specific azimuthal "quantum" number q are calculated as eigenvectors of the KL 
equation [201, (25)] 

R,{x, x')K^'^^ {x')dx' = A2_,if {x) (8) 



2I2; Jo 



with eigenvalues ^. The integral kernel is given by 



2tMx,x') = -c^a:' / ^^^^ + x^' - 2xx' cosOT^-^^'^de' + 2^x^[gi(x) + Giix') - ^2], (9) 







.2^ 1 

2t^Rq{x,x') = -c^x' j (x^ + .x^^ - 2a;x' cos6>0^^+^^/^ cos{q0')d0', q ^ 0, (10) 

( 2iL^_ 1 







where 



^1(2;) = c^-^-^^ f^^ dx"x" f^\x^ +x"^ -2xx"cos0")^^+^'>^^de", (11) 



Jo 

= (^/'^ci."."^r(."). (12) 

Since we do not assume that the expectation value of the tip-tilt component of the phase over the pupil vanishes, 
the terms proportional to Fried's variables C/3 and C/4 do not appear in the covariance and vanish in our analysis. 
Consequently, here and in [H, [5l|, tip-tilt modes are in the list of KL eigenfunctions, but not in Fried's list. 

B. Implementation 

In numerical practise, the KL equation is solved for the symmetrized variables, here marked with a tilde. 



k'f\x)^^^^K^^\xl (13) 

x' ~ 

Rq{x,x') = —Rg{x,x'), Rq{x,x')^Rq{x',x), (14) 



which turn (ISl) into 



xx'Rq{x,x')kl''\x')dx' = Xl^qkl^^x). (15) 



The benefits of working with an integral kernel that is symmetric under the exchange x ^ x' are 

• numerical stability and performance by use of linear algebra eigenvalue solvers for this type of symmetry. 

• immediate evidence that the eigenvalues are real- valued with orthogonal eigenvectors, 

Jo 



k\^\x)kl;!\x)dx = V: (16) 



xK'f\x)Kl^\x)dx = V- (17) 
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Further comments on the numerical treatment are given in Appendix IA II Some variants of the actual representation 
of the radial functions Kj,^^ exist. A simple and stable format is the finite-element (FEM) representation, in which 
Kp'^'' is a vector of values on (possibly equidistant) Xi, i = 1, 2, . . . , iV. In this case the matrix representation of the 
kernel of the KL equation (fT5|l is the N x N table of the gridded \/ xx' Rq{x , x') multiplied by any weights associated 
with the numerical integration. A power basis 

oo 

i^(«)(a;)=xi/2+l9l^fc,p,a;^- (18) 

3=0 

has the (small) disadvantage that — after insertion into (jlSp and out-projection of the kipq — a non-diagonal overlap 
matrix is left on the right hand side which leads to a slightly more complicated generalized eigenvalue problem. This 
is easily avoided by moving on to a Zernike basis of orthogonal Jacobi Polynomials. 

Figures of the basis functions F{r) with the largest eigenvalues, which are dominant and represent the speckles of 
largest size have been shown before f5ll, Fig. 5b]. The order of the rotation axis is determined by the eigenmodes 
show up in pairs that can be mutually converted by rotation around the pupil center by angles of tt/ (2q), as established 
by Mq{e). 



C. Wave number (Fourier) Representation 



The two-dimensional Fourier transform of the reduced KL eigenfunctions Kp"^^ Mq is 

Fj{ct) = / (fxK^p''\x)Mq{e)exp{2TTicT-x)= xdx d9K^p''\x)Mq{9) exp[2TTiax cos{e„ - 9)] 

Jx<l^ Jo Jo 

= K(^\a)Mq{9,), (19) 
where a = |cr| and 9a- define the spherical coordinates of the wave number, and 

iv:^')(cr) EE 27ri« / xKl'^^{x)Jq{2TTax)dx; X^'^ (x) = 27r(-i)« / aK^p^\a)Jq{2Tiax)da (20) 



is a Fourier pair. Not to insert factors of Ij, in the definition of the Fourier transform here is a judicious choice to 
ensure that the normalization p7p is the same in the x and in the a domain: 

aif(^)(a)i^(?)*(a)da = V- (21) 

If the Kp'^\x) are expanded in a series of Zernike polynomials (App. [Cjl. the Kp'^\a) are the equivalent series of 
Bessel Functions [13, (8)]0. 



III. TAYLOR MODEL 



A. Equation-of-Motion and Gradient Matrix 



The theme of this paper is how any two of the sets of coefficients aj are morphed if time is a added as a new 
parameter to the description. 

To ensure steady transformation in short time intervals, we will employ the Taylor model of lateral displacement 
into the direction of a velocity vector v, which is represented by [l^ 

^{r,t)=^{r-vt,0) (22) 

for time t, distance r to the pupil center and azimuth 9 along the pupil periphery [4l[. We make the scaling of the 
coefficients aj by the eigenvalues A explicit by writing the expansion as 

(^(r,^,i) =^/3^'"^(i)AL™if/'"^(x)M,„(0), x^21,r/D. (23) 

l,m 
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The fundamental scaling parameters (D/ro)^'^'^ and A can all be partially absorbed in these factors or basis functions. 



The Qj in ([T]) obey a statistics with variances {D / tq)^^'' \f , so the components of the random vector introduced 
here have all the same variance, (D/ro)^^"' . The Taylor model is basically a means to substitute the time derivative 
in the equation-of-motion (EOM) of (/s by a gradient, dt — > — v • Vr, 

a/3(™)(-J-\ 

dMv.t) = ' \mK["'\r)M^{e) = -V • Vr^(r,0). (24) 

This is the infinitesimal version of de-centering the basis functions pH [23l . [29| . The corresponding requirement in 
Fourier space is 

ip{v, a) — if{cr)5{v — a ■ v). (25) 

In polar coordinates in the pupil plane, 

9 = arctan ^; r = ^/x"^ + Y'^, (26) 
the two components of the gradient operator are 

dx - |^ae + |^a, = -— 9,+cos0a„ (27) 
aX oX r 

= ^de + ^dr^^'^de+s^nOdr. (28) 
oY oY r 

To tighten the notation, we assume that the velocity has no component in the F-direction of coordinates, so (|28p is 
not needed and ([M)) becomes 



E '^^^^\mK[^\r)MUe) = -vdx'p{r,0) = v (^de - cosOdr) Y^(5f\t)\,.^Kf {T)M^{e) (29) 

i,m j-n 

R^'^hi\hf:^ll!fl ^» / -n sin 61 sin ne* flH^.N^il^^^El!!^!! v fcos6lcosn6l , , 

l^Pj W ""TT^i nsin^cosn^ ''Z^'^j rr^ 77^ 1 cos0sinn0 ^^^^^ 



J," 



.y27r 1 sin ^ cos nf? ^ j ^ ' Qj. -v/27r 1 cos 6 sin n6' 

(n) ... >'j,nK)"'{r) ( cos(n + 1)9 - cos(n - 1)9 



j.n>0 



2ti \ sin(n + 1)6^ — sin(n — 1) 



M(, M.dKf{r) ^„ r cos(n + l)0 + cos(n-l)0 
2Z^^J ^ ^ 5r \/2:^ 1 sin(n + 1)6* + sin(n - 1)61 • 



This write-up is a composite of an upper line for the even and a lower line for the odd m. The upper line refers to the 
Mm cosine modes, including the radially symmetric m = modes, and the lower line refers to the Mm sine modes. 
The coupling is between azimuthal parameters that differ by one, m ^ n ± 1, similar to the selection rules of the 
electric dipole operator between hydrogenic states of the Schrodinger atom. 



B. Hybridization of KL Eigenmodes 

To isolate one coefficient, we multiply (PT|) by a general \k,sK\^\x)Mk{9), substitute r xD / {2'tx), and integrate 
over the pupil, exploiting the orthogonality relations (jl7p : 

1. 



j-,n>0 ^ 

A . n 



n— i ,s 



i,n>0 
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J2 Z?]"^ / ^xdxA,,„9.ifj")A,,.if(^)^5„+i,. 

J2 /3f ' [^^ xdxA,-„a,iff )a,,.4^)^<5|„_ii,,], s + 0. (32) 



j,n>0 

This is a system of linear homogeneous differential equations with a skew-symmetric, real-valued, sparse matrix 51, 



^|/3i^'=E"'^--'4"'' ^^ = ^' = -^, (33) 



where the symbols T and f mean transposition and Hermite conjugation, respectively. Computation of the matrix 
elements in Fourier space is proposed in Appendix [Bl Existence of a skew-symmetric representation is expected from 
the fact that the gradient operator changes sign with the parity of a polar vector. In the Cartesian X-Y coordinate 
system, it is coupled to the sign change of the derivative integrals after partial integration (assuming the "surface" 
integrals of Green's law vanish) , but it is less obvious for the set of integrals reduced to the radial coordinates r or 2: 
and assembled in fl. The analytic proof works with the derivative of Mercer's theorem of the covariance function in 
the KL kernel, 

-^R,{x,x') = 5] A2,^if(«)(x)4«)(x') (34) 

and is omitted here. It is mainly to retain this symmetry feature of the gradient matrix that I chose to split off the 
values of A in ([M]) instead of working with the aj expansion coefficients. 

Skew-symmetry implies that the eigenvalues of the transformation to principal axes with some orthogonal matrix 

^ = iiojii^-, e^t^i, (35) 

are pairs of complex-conjugate, purely imaginary numbers iioi (sl. [sl [27l. [35l l48j . This transformation of the basis of 
the /3 vector induces a transformation of the KL basis, 

^ = ^ (Z?]"))'^ A.^^xj-^M™ = (r/3r^)''cA,,„ifj™)M„. (36) 

Whereas the KL eigenf unctions have a rotational symmetry in the pupil plane inherited from M(Q\ each gradient 
eigenfunction shows some undulation along the v direction (here: the direction of X) at a spatial period of 'kD/{wi). 

From the mathematical point of view, the standard aim of this basis transformation is to decouple the EOM's of 
the expansion coefhcients (j33p . 

Y^;d^ = ^"'(^^^ ^^^^ 

which solves the time-dependent differential equation in terms of oscillatory fluctuations from initial values, 

Pi{t)^exp(^^ioJi?j Pm^o- (38) 

^ (s) 

where /3; is the l-th component of the matrix- vector product . These solutions at discrete angular frequencies 

CO = 2vuJi/D explain in conjunction with (|25p why the diagonalized eigen-modes have well-defined spatial frequencies 
along the v-direction. 

From a less formal point of view, the basis transformation is the definition of waves traveling in v direction under 
the conditions of 

• compatibility with the structure function, 

• enabling smooth dragging of the phase screen (compatibility with the Taylor hypothesis) by pairing of gradient 
eigenf unctions that are shifted relative to each other along w by a quarter of the spatial period. 
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FIG. 1: The solution (I38p of the differential equation places the components (here: m and I) of the vector p{t) of expansion 
coefficients on circular orbits with different velocities (bold arrows). To ensure that each component samples its full Gaussian 
distribution over large time intervals, a stochastic term needs to be added such that the simulated motion becomes chaotic as 
t increases (dashed circles with variable radii). 



This is to be compared with the alternative of starting from the structure function as a power density function in 
wave number space [H, IsJ. A key ingredient of the KL equation is its dependence on the finite support (here: the 
circular shape of the pupil) which samples the covariance 39]. The Gibb's oscillations which represent this cutoff are 
a major element of the Fourier decomposition of the KL functions. Working in real space with the gradient matrix 
might be interpreted as a deconvolution of these, followed by a superposition such that components on ridges of a 
given projection in the v-direction, i.e. aMiiOa) =const in our notation, are accumulated. 

The phase screen is generated by using these time-dependent coefficients as multipliers for the gradient eigenfunc- 
tions. The latter are static and have to be generated only once for a particular set of KL functions, 

T 

^\k,sK['\T)AUO) =Y,mFi{v), (39) 

I 

where Fi is the Z-th component of the matrix- vector product S^XjFj. The transformed initial coefhcients {£,* l3)\t=o are 
a linear superposition of Gaussian random variables and therefore Gaussian random variables themselves. Since ^ is 
orthogonal and since splitting off A in made the fi identically distributed (iid) , the independence is sustained by 
the transformation [36]. 

An accurate solution along these lines for some initial values does not actually produce a simple shift as a function 
of time; one reasons is that any finite basis set remains incomplete and the coupling to the fine-grained modes is 
missing — which will be discussed in Section ITVl 



exp 



IV. RESUMPTION OF STOCHASTIC TIME DEPENDENCE 

The use of the Taylor hypothesis in Section IfffI puts a tangent condition on the components (3i it) at each point in 
time, which lets them wander on circles in the complex plane once some randomly selected A|t=o have been chosen. 
For a finite set of basis functions, say 1 < j < in ([1]), we have essentially constructed some smooth embedding of the 
initial, randomized phase screen into an infinitely large phase screen, which turns out to be periodic in N variables. 
This determinism looks like an artifact if one aims at creation of realistic time series; actually, each component ought 
become independent of (uncorrelated to) its initial value on large time scales, which calls for some modification of the 
motion of the components as indicated in Fig. [TJ This would possibly be achieved by conversion of the deterministic 
differential equation ((24)) into some Ornstein-Uhlenbeck process by adding a time-dependent stochastic term 

We conclude with remarks on how the Taylor screen ansatz can be employed to dither the time-dependent expansion 
coefficients without searching for such an explicit (in the nomenclature of differential equations: inhomogeneous) term 
of the EOM. The differential form is not strictly equivalent to the translational symmetry (P^ : The value of 
(/?(r — vt,0) on the right hand side of ((22)) is generally — after times of the order D/v — peeking outside the pupil of 
size D\ it is not known precisely at i = because the N basis functions and associated expansion coefficients have 
only been endowed with the representation of the phase screen inside D. With the Taylor hypothesis we can virtually 
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reach out to r — vt if we move on to a size 



D > |r- vtl 



embedding it such that in the inner part consistently 



ipD{r,t) ^ ifoir - vi, 0) = v?£,(r - vi,0), 



(40) 



(41) 



where ipjy is a hnear adaptive function of (po ■ A set of higher order adaptive functions within a super-aperture D > D 
would contain the full information to implement (|22[) at finite times, and would be constructed with a basis set size 
N > N to account for the additional degrees of freedom to represent the speckles in the annulus outside D. 

For a quantitative model, one assumes that an excess N — N basis functions are discarded while synthesizing the 
phase screen. Working with the smaller basis set at the argument x = r/D, the rightmost N — N columns of the full 
matrix transformation from the arguments r/D of the host phase screen 



■if^«)(21,ri/Z?) K[''\2t,n/D) 



K^^\2t,r2/D) K^^\2t,r2lD) 



Kj^\2t,n/D) 

Ki''\21,r,/D) 



' ^ 00 

n(9) 
L 01 



v 



^ 10 
p(9) 
^ 11 



p(9) -r(9) 

^ Ofc Ifc 



(42) 



would have been discarded. So it is the product of these bottom rows of the F matrix by the K^^"* of higher order A:, to 
be multiplied by coefficients /3^, that has been ignored in the analysis of Section [1111 To re- insert their time-dependent 
randomness into the simulation, one can build the dot product of the missing rows and columns in (I42p . and multiply 
this vector with coefficients /3 that are randomly generated at each new image on the fly, consistent with the fact 
that they have unknown history, that they are not carried over between the time steps of the simulation, and that 
they have no associated eigenfrequencies in the gradient matrix. The elements of the components of T{D/D) depend 
on the ratio D/ D and scale these random contributions appropriately depending on how large the size D of the 
super-aperture must be chosen to comply with (|40p during the next time interval. 

In conclusion, feeding a stream of random kinks into the EOM of the expansion coefficients that otherwise move on 
circles does not introduce more parameters; it can be derived from a model of coupling the information just outside 
the pupil into the actual diameter once the effect of radial scaling on the basis functions has been calculated. 



V. SUMMARY 



The numerical synthesis of — in the Markov sense — independent phase screens by superposition of statistically 
independent basis functions over a telescope entrance pupil multiplied by a vector of Gaussian random numbers is a 
known concept. Connection of these still pictures to a movie needs some specification of the equation-of- motion of the 
random numbers, which I have bound to the Taylor "frozen" screen assumption that on short time scales an instance 
of a phase screen is merely shifted. 

If one external parameter — the wind velocity — is introduced, this suffices to formulate a first-order differential 
equation in time for the deterministic ("ballistic") motion. A decoupling of temporal and spatial structure is found 
by diagonalizing the skew-symmetric matrix of the gradient operator. This diagonalization needs to be done only 
once for each size of the basis set (each Taylor movie), and introduces new basis functions which undulate across the 
pupil in the wind direction on scales determined by the eigenvalues of the diagonalization. 

Randomness in this scheme of oscillating expansion coefficients (in the diagonalized coordinate system) is re- 
introduced by virtual embedding of the pupil into a larger pupil which regulates at which time scales and at which 
strength the fluctuations outside the pupil — hence unknown — mix higher-order fluctuations into the time-dependent 
expansion coefficients. 

The benefit of this technique is in the reduction of the problem of two-dimensional time-dependent phase screens 
to the time-dependence of a vector of a few expansion coefficients. This aims at more efficient implementation on a 
computer and is equivalent to pulling a much larger sky linearly across the circular telescope's aperture. No new physics 
is introduced; on the contrary, the technique is largely independent of modal decomposition and parametrization of 
the phase structure function. 
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APPENDIX A: KL INTEGRAL KERNEL 
1. Recursion with respect to azimuthal quantum number 

The kernel of (jS]) and P^. 

~c^— V^S,{x,x')K'f\x')dx' = \l^xK^^\x), (Al) 



21. JO 



contains functions Sq of the form 



Sq{x,x)^— — -2—/ de^/^{x^ ^x''^ -2xx' COS 9)''^+''^''^ COS qO. (A2) 
Auxiliary polar coordinates 

a; = 21a;UC0sa, x =2^xUS\n-a, sin(2Q) = 2;, (A3) 
render this into hypergeometric functions 

Sq = u'+''^-j'J de{i -zcoser+'<y^cosqe = u'+'< ^2ntlIl_2Mi ^ly ,F,(^'^~ + Iz^j 

. . ^^^^ 

in the argument range < z < 1- The quadratic transformations of 1, 15.3.19-21] are applicable or a re-writing as 
Legendre functions of [T, 15.4.10-11], which may be used to tune into whatever numerical library is available. 

2. Fourier Transform 

The Fourier representation of the KL equation (|A1|) is obtained by switching the {x,x') coordinates in Sq to the 
wave number domain introduced in Section III CI We summarize the transformation for q ^ 0. First, the quadratic 
transform [1,, 15.3.20] is applied to the hypergeometric function (|A4p 

(A5) 



Sq = u^^-'.f^^.t^I^-A (1) " (1 + 1/2 - 7/2, q + 1/2 _2^\ 

' V 2 <?! V2/ ^ ^ ' 2q + l ^l + zj 



Back to the a-coordinates with z — 2 sin a cos a and 1 + z ~ (sin a + cosa)^, 

Sq = u'+'^{smacosar'^'2ntlIl_lMl(^l + ,)-«+(i+7)/2 ^pf<l~ 1/2 - 7/2, q + 1/2 4sinacosa \ 

q\ \ 2q + 1 (sma + cosa)"' / 

_ (A6) 
its argument now matches the integral [22, 6.576.2] provided q > (1 + 7)/2 

^-2-7 j^^g-j^ afc) Jq(cos ak)dk 

- (cosasina) ^,^ , . ^^'^ ' f 7 /o /o^ ( ' ' " ' ^ I (A7) 

^ ^ 22+7(sina + cosa)29-i-7g!r(3/2 + 7/2) \^ 2g + 1 ' (sin a + cos a)2 ^ ^ ' 

The equivalent write-up 

2'+^ j^(^/^ + 7/2) r^-^-rj (sin ak)Jq (cos ak)dk 
1 (-1/2 - 7/2) Jq 

^zy (-1/2-7/2)^(1 + ^ /^g-l/2-7/2,g-H/2 , 4sinacosa ^ 

-I2J ql ^^M 29+1 '(sina + cosa)2 I ^^^^ 
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allows replacement of the hypergeometric function by the fc- integral in (jA5|) . 

Sg = M^+^(sinQ;cosa)^/^27r2^+^ / k-^-'' Jg{smak)J„{cosak)dk. (A9) 

1 (-1/2 -7/2) Jo 

We substitute (sin a cos a)'^/^ = (xa:')^/^/(21j;u) and insert Sq back into (jAl[) . 



1 .1. 



da;^a;x'it^+^7r2'^+^ w^^^l^ '^^^1 -^^"^^ (^0 / k-^'''Jq{smak)J„{cosak)dk = Xl„xKl'i^{x). (AlO) 
1 (-1/2 -7/2) 7o 



^'^ (21^)2 y„ ---^ - r(-l/2-7/2)^ 
The linear substitution k = An'txucj of the integration variable turns this into 



r-1, 

^''W^f+^J, - - r(-l-7/2) 



-c^ .^. / dx':Ex'7r22^i^^/^±^X^'?)(x')7r-i-^ H o-^-\Jq{2^xa)Jq{2^x' o)da ^ \1,jXKI^\x). (All) 



The momentum representation of the KL equation is obtained by multiplication with Jq{2Tia'x) on both sides and 
integration over x with ([^0]) . 

1^+7 (2^)' +7 r^%y-Ty^/2) ^^^'"'"-^^(^^^^O J.(2^^^)4^H^) = Ag,,4'')(aO, (A12) 

where 

r(3/2 + 7/2) 
^'^(27r)i+7r(-l/2-7/2) 

is a well-known factor [H (3)][3l, (25)] at 7 = 2/3. 



0.022656 (A13) 



3. Connection to Noll's Covariances 

If the radial functions are expanded in an orthonormal Zernike basis (sil . l33 | , 

?9 r^l-n ^J "n.Ti'Ji^an 



xi?,«,(x/l.)i?,^,,(x/l.)dx = -f^; (A14) 



K'f\x)^ J2 r„,p,,i-v/2(n + l)i?^(x/l,), (A15) 

(lAlip becomes 



1. 

n— g(mod2) 



vr r(3/2 + 7/2) f " ^^/^.^z ^ ^^3(71 + l)i?« (a;'/!^) / (7-2-7 j_^(27ra;CT) J,(27ra;'CT)da 

^ XlgxJ2rn',p,,VW+^K'i^/^^)- (A16) 



^21+713+7 r(-l/2- 7/2) io 



Integration over x' on the l.h.s. with the Fourier representation (|^n|) [tI. [qI. [TtI. [33| 

Vr„,p,,2«V2(«+l) xRUx/l.)Jq{27rax)dx ^ V T„,p,,z"l2 ^2(n + 1) ^"+1(2^^^^^) (A17) 

n n 

yields 

.-^ r(3/2 + 7/2) y2(;rTT)z"-ll r^^^±^^^^a'^-^J(2.xa)da 

^21+713+^ r(-l/2- 7/2)^ ^ "io 2w Jg(^^xa)«a 

= Xl,qxJ2r„,,p,q^2in' + l)Rl{x/l,). (A18) 
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Multiplication by R'^„{x/tx) on both sides and integration over x employing (jA14|) on the r.h.s. and integration over 
X with (|A17|) on the l.h.s. yields 

TT-'f r(3/2 + 7/2) ^ J„+i(27ral,) _2_^ J„"+i(2wl,) 
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The substitution 2TTaTLx — k of the integration variable generates 
r(3/2 + 7/2) 



/•OO 

Vr„.p.,V2(^rTT)r-^ / dfcJ„+i(fc)fc-4-^r"-V„..+i(fc) = , 

"'0 



(A20) 



r(-l/2-7/2)^'"^^^''^^^'^ --^^^^y-'- ■''■V2K + 1)' 

For n + n' > 1 + 7 > 0, the integral is [13, 6.574.2] 

where [1, ill 

r(4 + 7)r( "+"-^-^ ) 

Inn' — 7 ; \ 7 \ 7 \~ ■ (A22) 

24+7P I 2 + " -"+1+7 j P I 3 -|- "+" +1+7 J P / 2 + +1+7 j 

This is consistent with Noll's equation (25), the Kolmogorov covariance of a field with unit amplitude in terms of 
an orthogonal Zernikc basis {Zj}, 

< >= + ( —) / (A23) 

re- interpreted with (|A13p . 



This diverges from (|A2ip by a nominal factor 4/7r, 

• the factor 4 means that the Tn,p,q are amplitudes to be post-multiplied by D^, Fried's equation (23), whereas 
(IA23P contains an implicit factor of the squared radius which was silently dropped between Noll's (22) and (25). 

• the unit amplitude in Noll's Zernike basis — which is normalized with an additional weight I/tt in his (3) but 
not in my ([TT]) — is equivalent to an amplitude ^/tt here, and introduces a factor tt because the covariance is a 
quadratic form of the amplitudes. 

On the same basis, Dai's variable cq Jii, (2.10)] is consistent with the present calculation. Roddier's expression 
E{aj,aji) contains an additional factor 2^+^/7r « 1.0106, and the same small factor is present in Conan's version 
of (|A13p [13i, (3.20)] (U, (15)]. I am unable to trace these back — if applied it triggers a systematic increase of all 
eigenvalues ,j by one percent. 



4. Piston Term at g = 

The case of the modes Kp'^^ in ([8|) induces three special terms proportional to Gi and Q2 after insertion of ([9]). 
Modes on both sides of the eigenvalue equation are expanded in Zernike bases with (jA15|) . 

Isolation of a single r on the right hand side, ie, a matrix equation for the r, is achieved by multiplication with 
xK^\x) and integration over dx. A special case of (jA14[) . namely 



pi 



5 1^ 

xdxRlix/t^) = / xdxRl{x/t^)Rl{xlt^) = 
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reduces some integrals at the Gi and G2 to yield 



—CipSir Tpfi^n'Vn' + 1 / u^^^ du da zR'^{2u cos a) R'^^,{2u sin a) 2 Fi j 

n'=0(inod2) \ 

+<5Ti,oCi^87r Tpfi^n'Vn' + 1 / u'*+'''(iw(iazi?°, (2ucosq;) 2-^"! I 

n'=0(mod2) V 



-1/4 -7/4, 1/4 -7/4 



1/4 -7/4, 1/4 -7/4 



+Sn', oTpfl^oc^Sn J u^+'^du da zR^„{2u cos a) 2Fi\^ 7/4,1/4 7/4 |^2j 

TT 1 
— -7^2'5n,0'5n'.0''p.O,0 — o'''p,0.n , (A25) 

4 ' yyn+ 1 

The value of (fT2|) is a constant, the residual of Dai's first mode [H Table I]: 

^2 = c^8^ f ^ dxx f dx"x" f de{x^ + x"^ - 2xx" cosey^+'^^Z^ (A26) 

7r/4 (21j;)''+T Jo Jo 

which is transformed to the polar coordinates proposed in (jA3p . with Jacobian {21x)'^u, 



1 



1^ (.li. (.7r/4 (.l/(2cosct) (.7r/2 (.l/(2sinQ) 

dx / dx" = / da / udu + / da / wdit. (A27) 



(2ll.)2 

JO Jo Jo Jo Jtv/4 Jo 



r 16 

t/2 



//'7r/2 /'27r 
, udw / daw^ sin(2a)Mi+'^ / d6l[l - sin(2a) cos 6'l(i+'^)/2 (A28) 
0<usina<l/2 ^ ^ Jo 



TT J < u sin a < 1 /2 
< ucosa < 1/2 



'-^c, J u^^-du J daz{2n)2F, -1/4-7/4,-7/4+1/4 1^2^ (^29) 
64c. r J 2^r f -^/^ - - \A (A30) 







(5 + 7)(2cosa)5+T ^ M 1 

1 r(3 + 7) _ r(3/2 + 7/2) 

(5 + 7)21+^ r2(5/2 + 7/2) ''^r(-l/2-7/2) 



4cv7 '^^rnr- T.2fK /n /o^ 8cy / ^ ;^ /oo w 0.299953532617054c^ w 1.032421639. (A31) 



Since Rg — 1, within the top row 71 = of the coefficients matrix, the first two and the last two lines of (|A25|) cancel. 
Within the left column n' = of the coefficients matrix, the third line cancels the first and the fourth line cancels the 
second. This is the real-space equivalent to Noll's removal of the piston mode in Fourier space [33 |. which eliminates 
the Rq term from the q = block of the KL equation. 

APPENDIX B: GRADIENT MATRIX IN WAVE NUMBER SPACE 

The Fourier integral of the EOM p4)) transforms the gradient into a multiplication with the wave number —2mcT, 

E ^^^A,„if(™\^) = -V ■ (-2.0^^^/3(")(t)A,,„<)(^). (Bl) 

l.m j^n 

As before we assume the wind velocity v contains only a component into the X direction, 

V ■ (T = va cos9a, (B2) 

and HI]) yields 

E ^^^Ai,™K(™)(a)Af„,(0.) = 2™,^cTcos0.^/3f )WA,-„A'f )(a)M„(0,). (B3) 

l.m j,n 
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After insertion of ([6]), 



^d^^^ („,) r co.m0^ cos(n+l)0. + cos(n-l)^?. 

dt ' ^ '\J 2TT\sinm0^ D j^^^ ^ ^ ^ ^ ^ V 2^ \ sin(n + 1)0, + sin(n - 1)0, 

we project this onto a single component of the /3-vector by multipUcation with aXk.sK^f^^''* {a)Mk{0a) and integration 
over (Pa using the orthogonahty (pij) on the left hand side, 

_^rf^_^2 = ^^(")(^) 2 ^ r daa^\,,^Kf\a)\,^sK^:'*{a). (B5) 

21:r7riw at ^ ^ JenSs Jo ■' 

For a Zernike decomposition (|A17[) . the integrals Jj' +i{(t) Jk'+i{cr)da can be evaluated analytically [30, p. 50], 

o ^j.^ ^ ^k,s = y^/3-"''(t) ^ ((5n+l,a + ^Ti-l,a)(l ^ Ss+n,-l)^j.n>^k,s 

X 5](-l)(^'-'='+l^'-'='l)/V,,- „V(.?' + + (B6) 
j'fc' 

which recovers the elements of Noll's gradient matrix j^-, [33 |. 



APPENDIX C: ZERNIKE COVARIANCES 

A table of the dominant KL eigenvectors K^^ Mq is attached, diagonalizing a 80 x 80 matrix for each individual 
\q\. It contains the eigenvalue, after a colon the value of \q\, then the expansion in Zernike polynomials in Noll's 
nomenclature. For g = 0, the value in parentheses is the index of Z , for g ^ a comma-separated pair of first the 
index for the cosine-term, q > 0, then the index for the sine-term, g < 0. 

Terms with expansion coefficients down to a threshold of 10~^ are listed. The squared expansion coefficients sum 
to I/tt which implies that the square of the eigenvector integrated over the pupil is normalized to unity. 

0.3529054 : 1: +0 . 5639006*Z(2 ,3) -0 . 0180275*Z(8,7) +0 . 0010116*Z(16, 17) +0.0000047*Z(30,29) 
+0.0000011*Z(46,47) 

0.0187926 : 2: +0 . 5550339*Z(6 , 5) -0 . 1006468*Z(12 , 13) +0.0108338*Z(24,23) -0.0003738*Z(38,39) 
+0 . 0000221*Z (58 , 57) +0 . 0000013*Z (80 , 81 ) 

0.0187926 : 0: +0 . 5550339*Z(1) -0. 1006468*Z(4) +0.0108338*Z(11) -0 . 0003738*Z(22) +0.0000221*Z(37) 
+0.0000013*Z(56) 

0.0052209 : 3: +0 . 5415715*Z(10 ,9) -0. 1562750*Z(18, 19) +0 . 0241979*Z(32 ,31) -0 . 0016589*Z(48 ,49) 
+0.0001053*Z(70,69) 

0.0048842 : 1: +0 . 0175046*Z(2 ,3) +0 . 5382264*Z(8,7) -0 . 1661406*Z(16 , 17) +0.0266343*Z(30,29) 
-0 . 0019283*Z (46 , 47) +0 . 0001237*Z (68 , 67) 

0.0021547 : 4: +0 . 5274038*Z(14, 15) -0. 1965600*Z(26,25) +0.0387868*Z(40,41) -0 . 0038510*Z(60 ,59) 
+0.0003027*Z(82,83) -0 . 0000093*Z(110, 109) +0 . 0000018*Z(140 , 141) 

0.0016529 : 2: +0 . 0944936*Z(6 , 5) +0 . 4941255*Z(12 , 13) -0.2485822*Z(24,23) +0.0581153*Z(38,39) 
-0.0070710*Z(58,57) +0 . 0006213*Z(80 ,81) -0 . 0000294*Z(108, 107) +0 . 0000031*Z(138 , 139) 

0.0016529 : 0: +0 . 0944936*Z(1) +0.4941255*Z(4) -0.2485822*Z(11) +0 . 0581153*Z(22) -0.0070710*Z(37) 
+0.0006213*Z(56) -0 . 0000294*Z(79) +0.0000031*Z(106) 

0.0010850 : 5: +0 . 5136656*Z(20 ,21) -0 . 2270441*Z(34,33) +0.0534941*Z(50,51) -0 . 0068153*Z(72,71) 
+0.0006519*Z(96,97) -0 . 0000346*Z(126 , 125) +0 . 0000037*Z(158 , 159) 

0.0007763 : 3: +0 . 1423413*Z(10 ,9) +0.4455925*Z(18, 19) -0.3017149*Z(32,31) +0 . 0907617*Z(48 ,49) 
-0.0149565*Z(70,69) +0 . 0017093*Z(94,95) -0 . 0001269*Z(124, 123) +0 . 0000101*Z(156 , 157) 
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0.0007595 : 1: +0.0040371*Z(2,3) +0. 1508302*ZC8,7) +0 . 4388195*Z(16 , 17) -0.3064797*Z(30,29) 
+0.0938139*Z(46,47) -0.0157202*Z(68,67) +0.0018198*Z(92,93) -0.0001376*Z(122, 121) +0.0000109*Z(154,155) 

0.0006175 : 6: +0.5007017*Z(28,27) -0.2508042*Z(42,43) +0.0677902*Z(62,61) -0.0103942*Z(84,85) 
+0 . 001 1738*Z (112,111) -0 . 0000838*Z ( 142 , 143) +0 . 0000076*Z ( 178 , 177) 

0.0004293 : 4: +0 . 1743752*Z(14, 15) +0.3990822*Z(26,25) -0.3363293*Z(40,41) +0. 1220098*Z(60,59) 
-0.0249092*Z(82,83) +0.0034953*Z(110, 109) -0 . 0003387*Z(140 , 141) +0.0000288*Z(176, 175) 

0.0003857 : 2: -0.0317365*Z(6,5) -0.2140331*Z(12, 13) -0.3498546*Z(24,23) +0.3572055*Z(38,39) 
-0.1430459*Z(58,57) +0.0318299*Z(80,81) -0.0047925*Z(108,107) +0.0005028*Z(138, 139) -0.0000441*Z(174,173) 
+0.0000017*Z(212,213) 

0.0003857 : 0: -0 . 0317365*Z(1) -0.2140331*Z(4) -0.3498546*Z(11) +0.3572055*Z(22) -0 . 1430459*Z(37) 
+0.0318299*Z(56) -0.0047925*Z(79) +0.0005028*Z(106) -0 . 0000441*Z(137) +0 . 0000017*Z(172) 

0.0003822 : 7: +0.4885875*Z(36,35) -0.2697207*Z(52,53) +0.0814295*Z(74,73) -0.0144431*Z(98,99) 
+0 . 0018757*Z (128 , 127) -0 . 0001646*Z (160 , 161) +0 . 0000148*Z ( 198 , 197) 

0.0002627 : 5: -0. 1968506*Z(20,21) -0.3562481*Z(34,33) +0.3585383*Z(50,51) -0. 1508079*Z(72,71) 
+0.0363009*Z(96,97) -0.0060054*Z(126,125) +0.0007061*Z(158, 159) -0.0000683*Z(196,195) +0.0000037*Z(236,237) 

0.0002516 : 8: +0.4773002*Z(44,45) -0.2850218*Z(64,63) +0.0943101*Z(86,87) -0.0188400*Z(114,113) 
+0 . 0027559*Z (144 , 145) -0 . 0002835*Z (180 , 179) +0 . 0000270*Z (218 , 219) 

0.0002258 : 3: -0 . 0579989*Z(10 ,9) -0.2465801*Z(18, 19) -0 . 2674992*Z(32 ,31) +0.3804364*Z(48,49) 
-0.1874358*Z(70,69) +0.0513178*Z(94,95) -0 . 0094502*Z(124, 123) +0 . 0012400*Z(156, 157) -0 . 0001291*Z(194, 193) 
+0 . 0000090*Z (234 , 235) -0 . 0000012*Z (280 , 279) 

0.0002234 : 1: -0.0015503*Z(2,3) -0.0625126*Z(8,7) -0.2490043*Z(16, 17) -0.2604727*Z(30,29) 
+0.3814605*Z(46,47) -0. 1902056*Z(68,67) +0.0525678*Z(92,93) -0 . 0097567*Z(122, 121) +0.0012903*Z(154,155) 
-0 . 0001351*Z(192 , 191) +0 . 0000095*Z (232 , 233) -0 . 0000012*Z(278 , 277) 

0.0001735 : 9: +0.4667824*Z(54,55) -0.2975492*Z(76,75) +0. 1064060*Z(100, 101) -0 . 0234853*Z(130, 129) 
+0.0038066*Z(162,163) -0.0004455*Z(200,199) +0.0000458*Z(240,241) -0.0000022*Z(286,285) 

0.0001723 : 6: -0.2131220*Z(28,27) -0.3173321*Z(42,43) +0.3721298*Z(62,61) -0 . 1768017*Z(84,85) 
+0.0486138*Z(112,111) -0.0092134*Z(142,143) +0.0012625*Z(178, 177) -0.0001388*Z(216,217) +0.0000104*Z(260,259) 

-0.0000013*Z(306,307) 

0.0001445 : 4: -0.0803076*Z(14, 15) -0.2618432*Z(26,25) -0. 1948668*Z(40,41) +0.3859041*Z(60,59) 
-0.2255185*Z(82,83) +0.0727395*Z(110,109) -0 . 0157444*Z(140, 141) +0.0024554*Z(176,175) -0.0002991*Z(214,215) 
+0 . 0000271*Z (258 , 257) -0 . 0000028+Z (304 , 305) 

0.0001363 : 2: -0 . 0144405*Z(6 , 5) -0 . 1076272*Z(12 , 13) -0.2659802*Z(24,23) -0. 1497605*Z(38,39) 
+0.3840616*Z(58,57) -0.2422123*Z(80,81) +0.0825518*Z(108,107) -0 . 0187019*Z(138, 139) +0.0030449*Z(174, 173) 
-0 . 0003847*Z (212 , 213) +0 . 0000366*Z (256 , 255) -0 . 0000036*Z (302 , 303) 

0.0001363 : 0: -0.0144405*Z(1) -0. 1076272*Z(4) -0 . 2659802*Z(11) -0. 1497605*Z(22) +0 . 3840616*Z(37) 
-0.2422123*Z(56) +0.0825518*Z(79) -0 . 0187019*Z(106) +0 . 0030449*Z(137) -0 . 0003847*Z(172) +0 . 0000366*Z(211) 
-0.0000036*Z(254) 

0.0001242 : 10: +0.4569675*Z(66,65) -0.3079023*Z(88,89) +0 . 1177299*Z(116, 115) -0.0282996*Z(146,147) 
+0.0050164*Z(182,181) -0.0006545*Z(220,221) +0.0000729*Z(264,263) -0.0000048*Z(310,311) 

0.0001190 : 7: -0.2251569*Z(36,35) -0.2821207*Z(52,53) +0 . 3795853*Z(74,73) -0. 1999679*Z(98,99) 
+0.0614426*Z(128,127) -0.0130623*Z(160,161) +0 . 0020318*Z(198, 197) -0.0002511*Z(238,239) +0.0000229*Z(284,283) 
-0.0000025*Z(332,333) 

0.0000983 : 5: -0.0988260*Z(20,21) -0.2668319*Z(34,33) -0. 1319724*Z(50,51) +0.3797889*Z(72,71) 
-0.2570938*Z(96,97) +0.0949715*Z(126,125) -0 . 0235100*Z(158, 159) +0.0042219*Z(196,195) -0.0005896*Z(236,237) 
+0 . 0000636*Z (282 ,281) -0 . 0000065*Z (330 , 331) 
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0.0000917 : 11: +0 . 4477892*Z(78 ,77) -0 . 3165203*Z(102 , 103) +0 . 1283147*Z(132 , 131) -0.0332196*Z(164, 165) 
+0.0063723*Z(202,201) -0.0009130*Z(242,243) +0 . 0001100*Z(288 ,287) -0.0000089*Z(336,337) +0 . 0000014*Z(390 ,389) 

0.0000901 : 3: -0.0293974*Z(10,9) -0. 1411969*Z(18, 19) -0.2594816*Z(32,31) -0.0576767*Z(48,49) 
+0.3648102*Z(70,69) -0.2824450*Z(94,95) +0. il41543*Z(124,123) -0 . 0303990*Z(156, 157) +0.0058362*Z(194,193) 
-0 . 0008646*Z (234 , 235) +0 . 0000997*Z (280 , 279) -0 . 0000103*Z (328 , 329) 

0.0000895 : 1: -0 . 0007599*Z(2 ,3) -0 . 0319163*Z(8,7) -0 . 1439972*Z(16 , 17) -0.2582430*Z(30,29) 
-0.0521890*Z(46,47) +0.3632810*Z(68,67) -0.2842444*Z(92,93) +0 . 1156663*Z(122, 121) -0.0309709*Z(154, 155) 
+0.0059756*Z(192,191) -0.0008892*Z(232,233) +0.0001031*Z(278,277) -0.0000107*Z(326,327) 

0.0000855 : 8: -0.2341860*Z(44,45) -0.2502655*Z(64,63) +0.3825946*Z(86,87) -0.2204388*Z(114,113) 
+0.0744793*Z(144,145) -0.0174801*Z(180,179) +0.0030287*Z(218,219) -0.0004158*Z(262,261) +0.0000439*Z(308,309) 
-0.0000047*Z(360,359) 

0.0000700 : 6: -0. 1141626*Z(28,27) -0.2655091*Z(42,43) -0.0779168*Z(62,61) +0.3661379*Z(84,85) 
-0.2825223*Z(112,111) +0. 1171858*Z(142, 143) -0 . 0325258*Z(178, 177) +0.0065797*Z(216,217) -0.0010345*Z(260,259) 
+0 . 0001279*Z (306 , 307) -0 . 0000140*Z (358 , 357) 

0.0000694 : 12: +0.4391862*Z(90,91) -0 . 3237332*Z(118, 117) +0 . 1382027*Z(148 , 149) -0 . 0381955*Z(184, 183) 
+0.0078604*Z(222,223) -0 . 0012226*Z(266 , 265) +0.0001586*Z(312,313) -0 . 0000148*Z(364,363) +0.0000020*Z(418,419) 

0.0000634 : 9: -0.2410183*Z(54,55) -0.2213988*Z(76,75) +0 . 3823434*Z(100, 101) -0.2384139*Z(130,129) 
+0.0874939*Z(162,163) -0.0223892*Z(200,199) +0.0042602*Z(240,241) -0.0006422*Z(286,285) +0.0000761*Z(334,335) 
-0.0000084*Z(388,387) 

0.0000629 : 4: -0 . 0439206*Z(14, 15) -0. 1651472*Z(26,25) -0.2407349*Z(40,41) +0 . 0166192*Z(60 ,59) 
+0.3327821*Z(82,83) -0.3113005*Z(110, 109) +0 . 1454590*Z(140, 141) -0.0443771*Z(176, 175) +0.0097695*Z(214,215) 
-0.0016568*Z(258,257) +0.0002217*Z(304,305) -0.0000253*Z(356,355) +0.0000020*Z(410,411) 

0.0000607 : 2: +0.0077992*Z(6,5) +0.0613337*Z(12, 13) +0. 1781334*Z(24,23) +0.2255147*Z(38,39) 

-0.0478290*Z(58,57) -0 . 3169309*Z(80 ,81) +0.3197894*Z(108, 107) -0 . 1560387*Z(138 , 139) +0 . 0492359*Z(174, 173) 
-0.0111660*Z(212,213) +0 . 0019449*Z(256 , 255) -0.0002674*Z(302,303) +0 . 0000311*Z(354,353) -0.0000026*Z(408,409) 

0.0000607 : 0: +0.0077992*Z(1) +0.0613337*Z(4) +0 . 1781334*Z(11) +0.2255147*Z(22) -0 . 0478290*Z(37) 
-0.3169309*Z(56) +0.3197894*Z(79) -0 . 1560387*Z(106) +0 . 0492359*Z(137) -0 . 0111660*Z(172) +0 . 0019449*Z(211) 
-0.0002674*Z(254) +0 . 0000311*Z(301) -0 . 0000026*Z(352) 

0.0000536 : 13: +0.4311032*Z(104, 105) -0 . 3297940*Z(134, 133) +0. 1474395*Z(166, 167) -0.0431883*Z(204,203) 
+0.0094668*Z(244,245) -0 . 0015839*Z(290, 289) +0.0002200*Z(338,339) -0 . 0000230*Z(392 , 391) +0.0000029*Z(448,449) 

0.0000516 : 7: -0. 1269134*Z(36,35) -0.2602434*Z(52,53) -0.0315944*Z(74,73) +0.3476608*Z(98,99) 
-0.3023922*Z(128,127) +0. 1387961*Z(160, 161) -0 . 0425527*Z(198, 197) +0.0095407*Z(238,239) -0.0016640*Z(284,283) 
+0 . 0002307*Z (332 , 333) -0 . 0000274*Z (386 , 385) +0 . 0000023*Z (442 , 443) 

0.0000482 : 10: -0 . 2462062*Z(66 ,65) -0 . 1951759*Z (88 , 89) +0 . 3796826*Z(116 , 115) -0 . 2541152*Z(146 , 147) 
+0.1003177*Z(182,181) -0.0277123*Z(220,221) +0.0057268*Z(264,263) -0.0009388*Z(310,311) +0.0001226*Z(362,361) 
-0.0000143*Z(416,417) 

0.0000458 : 5: +0.0572248*Z(20,21) +0. 1816940*Z(34,33) +0.2158725*Z(50,51) -0.0753939*Z(72,71) 
-0.2938276*Z(96,97) +0.3299769*Z(126, 125) -0 . 1751494*Z(158 , 159) +0 . 0600933*Z(196 , 195) -0.0148722*Z(236,237) 
+0.0028343*Z(282,281) -0 . 0004291*Z(330,331) +0.0000545*Z(384,383) -0 . 0000053*Z(440,441) 

0.0000433 : 3: +0.0170071*Z(10,9) +0.0872520*Z(18,19) +0. 1958042*Z(32,31) +0. 1813467*Z(48,49) 
-0.1207320*Z(70,69) -0.2590231*Z(94,95) +0.3385998*Z(124,123) -0 . 1938146*Z(156, 157) +0.0703626*Z(194,193) 
-0.0182822*Z(234,235) +0.0036395*Z(280,279) -0.0005748*Z(328,329) +0.0000755*Z(382,381) -0 . 0000078*Z(438,439) 

0.0000431 : 1: +0 . 0004302*Z(2 ,3) +0 . 0185360*Z(8,7) +0 . 0894311*Z(16 , 17) +0. 1964790*Z(30,29) 
+0.1784645*Z(46,47) -0 . 1238988*Z(68,67) -0.2561964*Z(92,93) +0 . 3390955*Z(122, 121) -0. 1952323*Z(154, 155) 
+0.0711828*Z(192,191) -0.0185635*Z(232,233) +0.0037077*Z(278,277) -0.0005874*Z(326,327) +0.0000774*Z(380,379) 
-0.0000080*Z(436,437) 

0.0000422 : 14: +0.4234909*Z(120, 119) -0 . 3349000*Z(150 , 151) +0. 1560712*Z(186, 185) -0.0481676*Z(224,225) 
+0.0111783*Z(268,267) -0.0019968*Z(314,315) +0.0002955*Z(366,365) -0.0000339*Z(420,421) +0.0000042*Z(480,479) 
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0.0000391 : 8: -0 . 1375743*Z(44,45) -0.2525093*Z(64,63) +0 . 0080638*Z(86,87) +0 . 3262011*Z(114, 113) 
-0.3173578*Z(144,145) +0. 1594031*Z(180, 179) -0.0533561*Z(218,219) +0.0130941*Z(262,261) -0.0025033*Z(308,309) 
+0 . 0003830*Z (360 , 359) -0 . 0000495*Z (414 , 415) +0 . 0000049*Z (474 , 473) 

0.0000375 : 11: -0.2501393*Z(78,77) -0. 1712883*Z(102, 103) +0 . 3752347*Z(132, 131) -0.2677650*Z(164,165) 
+0.1128295*Z(202,201) -0.0333760*Z(242,243) +0.0074239*Z(288,287) -0.0013126*Z(336,337) +0.0001867*Z(390,389) 
-0 . 0000232*Z (446 , 447) +0 . 0000020*Z (508 , 507) 

0.0000344 : 6: +0.0691358*Z(28,27) +0 . 1926538*Z(42 ,43) +0. 1883543*Z(62,61) -0. 1211123*Z(84,85) 
-0.2517239*Z(112,111) +0.3399261*Z(142,143) -0 . 2023814*Z(178, 177) +0.0770075*Z(216,217) -0.0211096*Z(260,259) 
+0.0044561*Z(306,307) -0.0007503*Z(358,357) +0.0001051*Z(412,413) -0.0000118*Z(472,471) +0.0000014*Z(534,535) 

0.0000337 : 15: +0.4163054*Z(136,135) -0 . 3392075*Z(168, 169) +0. 1641421*Z(206,205) -0.0531094*Z(246,247) 
+0.0129820*Z(292,291) -0.0024606*Z(340,341) +0.0003862*Z(394,393) -0.0000479*Z(450,451) +0.0000060*Z(512,511) 

0.0000321 : 4: +0.0267586*Z(14, 15) +0. 1087334*Z(26,25) +0.2017551*Z(40,41) +0. 1340015*Z(60,59) 
-0.1709979*Z(82,83) -0. 1972564*Z(110, 109) +0 . 3432679*Z(140, 141) -0.2272400*Z(176,175) +0.0933061*Z(214,215) 
-0.0272720*Z(258,257) +0.0060951*Z(304,305) -0.0010833*Z(356,355) +0.0001592*Z(410,411) -0.0000191*Z(470,469) 
+0.0000022*Z(532,533) 

0.0000313 : 2: +0 . 0047026*Z(6 , 5) +0 . 0382423*Z(12 , 13) +0. 1210695*Z(24,23) +0.2006325*Z(38,39) 
+0.1125597*Z(58,57) -0. 1854229*Z(80,81) -0. 1762032*Z(108, 107) +0 . 3424840*Z(138, 139) -0.2359329*Z(174, 173) 
+0.0996008*Z(212,213) -0.0297915*Z(256,255) +0.0067958*Z(302,303) -0.0012312*Z(354,353) +0 . 0001841*Z(408,409) 
-0 . 0000225*Z (468 , 467) +0 . 0000026*Z (530 , 531) 

0.0000313 : 0: +0 . 0047026*Z(1) +0 . 0382423*Z(4) +0. 1210695*Z(11) +0.2006325*Z(22) +0. 1125597*Z(37) 
-0.1854229*Z(56) -0. 1762032*Z(79) +0.3424840*Z(106) -0 . 2359329*Z(137) +0 . 0996008*Z(172) -0.0297915*Z(211) 
+0.0067958*Z(254) -0 . 0012312*Z(301) +0 . 0001841*Z(352) -0 . 0000225*Z(407) +0 . 0000026*Z(466) 

0.0000304 : 9: +0. 1465398*Z(54,55) +0.2432536*Z(76,75) -0 . 0420131*Z(100, 101) -0.3030268*Z(130,129) 

+0.3280619*Z(162,163) -0 . 1787481*Z(200, 199) +0.0647183*Z(240,241) -0 . 0172125*Z(286 , 285) +0.0035720*Z(334,335) 
-0.0005961*Z(388,387) +0 . 0000836*Z(444,445) -0.0000094*Z(506,505) +0 . 0000012*Z(570 , 571) 

0.0000297 : 12: -0.2530992*Z(90,91) -0. 1494654*Z(118, 117) +0 . 3694619*Z(148, 149) -0.2795742*Z(184, 183) 
+0.1249441*Z(222,223) -0.0393127*Z(266,265) +0.0093431*Z(312,313) -0.0017691*Z(364,363) +0 . 0002714*Z(418,419) 
-0 . 0000360*Z (478 , 477) +0 . 0000035*Z (540 , 541) 

0.0000273 : 16: +0 . 4095083*Z(152 , 153) -0 . 3428415*Z(188 , 187) +0. 1716943*Z(226,227) -0.0579956*Z(270,269) 
+0.0148663*Z(316,317) -0 . 0029741*Z(368,367) +0.0004929*Z(422,423) -0 . 0000656*Z(482 ,481) +0.0000085*Z(544,545) 

0.0000265 : 7: +0.0797010*Z(36,35) +0. 1994106*Z(52,53) +0. 1601663*Z(74,73) -0. 1560338*Z(98,99) 
-0.2088996*Z(128,127) +0.3425979*Z(160,161) -0 . 2266633*Z(198, 197) +0.0946223*Z(238,239) -0.0284013*Z(284,283) 
+0.0065649*Z(332,333) -0.0012137*Z(386,385) +0.0001862*Z(442,443) -0.0000235*Z(504,503) +0.0000028*Z(568,569) 

0.0000245 : 5: +0 . 0363113*Z (20 , 21) +0. 1259883*Z(34,33) +0. 1996065*Z(50,51) +0 . 0877660*Z(72 ,71) 
-0.2031331*Z(96,97) -0. 1363122*Z(126, 125) +0 . 3365674*Z(158, 159) -0.2553567*Z(196, 195) +0. 1170904*Z(236,237) 
-0.0379733*Z(282,281) +0.0093975*Z(330,331) -0.0018516*Z(384,383) +0.0003010*Z(440,441) -0.0000405*Z(502,501) 
+0 . 0000050*Z (566 , 567) 

0.0000240 : 10: +0. 1541205*Z(66,65) +0 . 2330979*Z(88 ,89) -0 . 0710791*Z(116 , 115) -0 . 2790161*Z(146 , 147) 
+0.3351005*Z(182,181) -0. 1966763*Z(220,221) +0.0764443*Z(264,263) -0.0218564*Z(310,311) +0.0048838*Z(362,361) 
-0.0008808*Z(416,417) +0.0001332*Z(476,475) -0.0000165*Z(538,539) +0.0000021*Z(606,605) 

0.0000239 : 13: -0.2552942*Z(104,105) -0 . 1294720*Z(134, 133) +0.3627115*Z(166,167) -0.2897374*Z(204,203) 
+0.1366035*Z(244,245) -0.0454611*Z(290,289) +0.0114734*Z(338,339) -0.0023126*Z(392,391) +0 . 0003798*Z(448,449) 
-0 . 0000535*Z (510 , 509) +0 . 0000059+Z (574 , 575) 

0.0000235 : 3: +0.0107543*Z(10,9) +0.0575014*Z(18, 19) +0. 1434604*Z(32,31) +0. 1897436*Z(48,49) 
+0.0490395*Z(70,69) -0.2179204*Z(94,95) -0.0977964*Z(124, 123) +0 . 3285876*Z(156, 157) -0.2685761*Z(194,193) 
+0.1292745*Z(234,235) -0.0436019*Z(280,279) +0.0111673*Z(328,329) -0.0022712*Z(382,381) +0 . 0003801*Z(438,439) 
-0 . 0000528*Z (500 , 499) +0 . 0000066+Z (564 , 565) 

0.0000234 : 1: +0.0002680*Z(2,3) +0.0117488*Z(8,7) +0 .0591127*Z(16, 17) +0. 1445874*Z(30,29) 
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+0.1887365*Z(46,47) +0.0460816*Z(68,67) -0.2187900*Z(92,93) -0 . 0948401*Z(122 , 121) +0.3278423*Z(154, 155) 
-0.2695379*Z(192,191) +0. 1302224*Z(232,233) -0.0440532*Z(278,277) +0.0113124*Z(326,327) -0 . 0023063*Z(380 ,379) 
+0 . 0003868*Z (436 , 437) -0 . 0000538*Z (498 , 497) +0 . 0000067*Z (562 , 563) 

0.0000224 : 17: +0.4030654*Z(170, 171) -0.3459037*Z(208,207) +0. 1787669*Z(248,249) -0.0628120*Z(294,293) 
+0.0168202*Z(342,343) -0 . 0035357*Z(396 ,395) +0.0006164*Z(452,453) -0 . 0000873*Z(514, 513) +0.0000116*Z(578,579) 

0.0000209 : 8: +0 . 0890428*Z(44,45) +0 . 2029997*Z(64,63) +0. 1324492*Z(86,87) -0. 1821037*Z(114, 113) 
-0.1668971*Z(144,145) +0 . 3393270*Z(180, 179) -0.2477570*Z(218,219) +0 . 1125023*Z(262 , 261) -0.0366365*Z(308,309) 
+0.0091872*Z(360,359) -0 . 0018463*Z(414,415) +0.0003077*Z(474,473) -0 . 0000427*Z(536 , 537) +0.0000054*Z(604,603) 

0.0000194 : 14: -0 . 2568812*Z(120, 119) -0 . 1111043*Z(150 , 151) +0.3552477*Z(186, 185) -0.2984316*Z(224,225) 
+0.1477708*Z(268,267) -0 . 0517665*Z(314,315) +0.0138018*Z(366,365) -0 . 0029460*Z(420 ,421) +0.0005147*Z(480,479) 
-0 . 0000769*Z (542 , 543) +0 . 0000093*Z (610 , 609) -0 . 0000013*Z (680 , 681) 

0.0000193 : 11: +0. 1605609*Z(78,77) +0 . 2224571*Z(102 , 103) -0 . 0959649*Z(132 , 131) -0 . 2547779*Z(164, 165) 
+0.3390096*Z(202,201) -0 . 2131084*Z(242 , 243) +0.0883640*Z(288,287) -0 . 0269784*Z(336 ,337) +0.0064472*Z(390,389) 
-0.0012473*Z(446,447) +0 . 0002022*Z(508, 507) -0.0000272*Z(572,573) +0 . 0000035*Z(642 , 641) 

0.0000191 : 6: +0 . 0453401*Z(28 ,27) +0. 1395506*Z(42,43) +0. 1919744*Z(62,61) +0 . 0447795*Z(84,85) 
-0.2210912*Z(112,111) -0.0789246*Z(142,143) +0 . 3210884*Z(178, 177) -0 . 2777758*Z(216 , 217) +0. 1408688*Z(260,259) 
-0.0501513*Z(306,307) +0 . 0135944*Z(358,357) -0.0029353*Z(412,413) +0 . 0005225*Z(472,471) -0.0000776*Z(534,535) 
+0.0000102*Z(602,601) 

0.0000186 : 18: +0 . 3969465*Z(190, 189) -0.3484772*Z(228,229) +0. 1853962*Z(272,271) -0.0675481*Z(318,319) 
+0.0188337*Z(370,369) -0.0041438*Z(424,425) +0.0007572*Z(484,483) -0.0001134*Z(546,547) +0.0000157*Z(614,613) 
-0 . 0000013*Z (684 , 685) 

0.0000182 : 4: +0 . 0175669*Z(14, 15) +0 . 0749570*Z(26 ,25) +0. 1578453*Z(40,41) +0. 1697010*Z(60,59) 
-0.0074973*Z(82,83) -0 . 2285922*Z(110 , 109) -0 . 0261598*Z(140 , 141) +0.3019316*Z(176, 175) -0.2923630*Z(214,215) 
+0.1586920*Z(258,257) -0 . 0595914*Z(304,305) +0.0169141*Z(356,355) -0 . 0038087*Z(410,411) +0.0007046*Z(470,469) 
-0 . 0001088*Z (532 , 533) +0 . 0000147*Z (600 , 599) -0 . 0000015*Z (670 , 671) 

0.0000178 : 2: -0 . 0030615*Z(6 , 5) -0 . 0254743*Z(12 , 13) -0.0850827*Z(24,23) -0. 1618826*Z(38,39) 
-0.1591569*Z(58,57) +0 . 0257668*Z(80 ,81) +0 . 2284503*Z(108, 107) +0 . 0073571*Z(138 , 139) -0 . 2934941*Z(174, 173) 
+0.2969140*Z(212,213) -0 . 1652728*Z(256 , 255) +0.0632820*Z(302,303) -0 . 0182638*Z(354,353) +0.0041754*Z(408,409) 
-0.0007833*Z(468,467) +0 . 0001227*Z(530, 531) -0.0000167*Z(598,597) +0 . 0000018*Z(668,669) 

0.0000178 : 0: -0 . 0030615*Z(1) -0 . 0254743*Z(4) -0.0850827*Z(11) -0. 1618826*Z(22) -0. 1591569*Z(37) 
+0.0257668*Z(56) +0 . 2284503*Z(79) +0.0073571*Z(106) -0.2934941*Z(137) +0.2969140*Z(172) -0. 1652728*Z(211) 
+0.0632820*Z(254) -0.0182638*Z(301) +0.0041754*Z(352) -0.0007833*Z(407) +0.0001227*Z(466) -0.0000167*Z(529) 
+0.0000018*Z(596) 



APPENDIX D: CORRIGENDA 
1. Wang and Markey 

Section IV in [5l| suffers from an inconsistent normalization of integrals over the pupil. This becomes first apparent 
in (29), where a division by ttD^ is missing on the right hand side. (Equation numbers here in Appendix ID] refer to 
the original papers [5l| and The physical units of W in (2) and of G in (28) are unity; the unit of (j) is radians, 

which means with (26), the j3 also carry the unit radians; there is a divisor missing with units of square meters to 
compensate the units of the differential dPr' in (29). There are two philosophies to correct this. The smaller amount 
of editing is needed if (28) is re-normalized such that the factor nD^ disappears — which resumes Fried's notation. 
Alternatively, the list of changes consistent with keeping (28) is: 

• Divide the right hand side of (29) and the left hand side of (30) through nD^. 

• Divide the middle term and the right hand side of (33) through ttD^ . 

• Replace the factor r in front of the integral in (34) by r' . 

• Divide the left hand side of (35) through tiD^ . 
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• Multiply the right hand side of (37) by tt. This demonstrates that the eigenvalues in the ^ column of TABLE 
III are compatible with those in [20, TABLE I]. 

A second problem seems to have been induced by the mixed use of real-valued azimuthal functions in (9) and complex- 
valued azimuthal functions in (32). To keep (32), (43) is to be changed to 



"\l{x)Kl,{x)xdx=\]l'^^ 1^-/; (Dl) 

, ifpT^p'. 



2. Fried 

Typographic corrections to [13] are: 

• /„ in (9) ought read /*. 

• (in' in (16) ought read /?*,. 

. (/3*/3;,) in (17) ought read (/3*/3„.). 

• Q2{x) in (27c) ought simply read Q2. Qa{x) in (27e) ought simply read t/4. 

• A minus sign is missing in front of Qq in (28). 

• The factor B{x') should read yV(a;') in the integral in the first line of (34). 

• In the first line of (36a), the term ~2xx in the square root ought read ^2xx' . 

• In (36b), a factor x is missing in front of Qy,{x'), and exp(i0) ought read exp(i0'). 

• In (36c), a minus sign is missing in front of the integral, and exp(jg0) ought read ex.]i{iq0')d9' . 

• The variables i3„ ought be squared in (37a) and (37b). 

• {\itD^)^i in (40) ought read (^TrD^)-!. 

3. Noll 

Typographic clarifications to [s^j are 

• The name Wiener is misspelled after (18). 

• [T{P + l)l2f in the denominator of (20) is to be interpreted as [r{(P + l)/2}]2. 

• the term ~2n in the exponent of the second line of (25) ought read —2m. 

• T[{P + 1/2] in the denominator of the unnumbered equation on page 211 is to be interpreted as T\{P + l)/2]. 
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